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Resilience of hybrid optical angular momentum qubits to turbulence 
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Recent schemes to encode quantum information into the total angular momentum of 
light, defining rotation-invariant hybrid qubits composed of the polarization and orbital 
angular momentum degrees of freedom, present interesting applications for quantum 
information technology. However, there remains the question as to how detrimental 
effects such as random spatial perturbations affect these encodings. Here, we demon¬ 
strate that alignment-free quantum communication through a turbulent channel based 
on hybrid qubits can be achieved with unit transmission fidelity. In our experiment, 
alignment-free qubits are produced with g-plates and sent through a homemade tur¬ 
bulence chamber. The decoding procedure, also realized with g-plates, relies on both 
degrees of freedom and renders an intrinsic error-filtering mechanism that maps errors 
into losses. 


Successful transmission of information from one location to another depends on the physical system carrying 
the information and how robust it is under the unavoidable influence of the environment. Among all the physical 
candidates for carriers of quantum information, light represents the most reliable one for long-distance communication, 
due to its low interaction with the environment and its high speed. The most popular degree of freedom of light 
exploited in quantum technologies is the polarization of single photons. However, polarization encoding requires 
a shared reference frame between sender and receiver [T], which represents a limitation, specially if one considers 
scenarios where the relative reference frame changes, as for example in satellite communication. This obstacle can be 
overcome with “hybrid” rotational invariant single photons states where qubits are encoded into both the polarization 
and orbital angular momentum (0AM) degrees of freedom [2]-[4]. Recent progress on the control of the orbital 0AM of 
light has proven it also as a good candidate for storage, transmission, and manipulation of quantum information with 
promising applications [SHin]. Meanwhile, quantum technologies with hybrid systems develop rapidly unveiling new 
applications in metrology m , and more recently in quantum key distribution through free space [12]. However, spatial 
properties of light such as 0AM are severely affected by spatial perturbations such as the optical turbulence present 
in the atmosphere [T3H22] . 0AM optical systems encoding classical information have recently started to be tested in 
real scenarios over distances of the order of kilometers with positive results and perpectives for the transmission of 
quantum information [23] . 

Since hybrid qubits exploit 0AM properties, they are also susceptible to atmospheric turbulence. This raises the 
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question of whether the use of this additional degree of freedom does indeed provide a practical advantage. Here 
we show that, in the weak turbulence regime, unity transmission fidelity of rotational invariant hybrid qubits can 
in principle be achieved. This robustness is due to both the logical encoding and the decoding procedure, which 
uses a “q-plate” device IS [241 and a projection onto the fundamental gaussian mode with a single-mode fiber. A 
preliminary study was reported in Ref. [4], where it was predicted that spatial perturbations mirror symmetric in 
0AM space should result only in losses, rather than qubit state errors. One can gain intuition of this intrinsic error 
correction procedure as follows: Hybrid qubits employ polarization and 0AM eigenstates |/), | —/). The wavefunctions 
of these 0AM states are equivalent up to a mirror reflection, so a spatial perturbation with this symmetry affects both 
components of 0AM equally. In the decoding procedure, the projection onto a zero-order 0AM mode transfers the 
encoded information to the polarization degree of freedom in a robust fashion: The effects of the 0AM perturbation 
would factor as they affected both 0AM states equally. 

Nevertheless, a turbulent atmosphere perturbs the spatial structure of a beam by tilting, skewing, displacing or a 
combination of all of these in a stochastic way. In general, such effects do not respect any symmetry. However, in 
the weak turbulence regime the perturbation appears as a multiplicative spatially-dependent phase term [T4j . This is 
enough to guarantee mirror symmetry, so that our decoding method selects only the homogeneous contribution 0AM 
terms, resulting in high transmission fidelity and mapping all the inhomogeneous contributions into losses. 

Here, we quantify these losses and perform an experimental demonstration of the fidelity resilience under atmo¬ 
spheric turbulence produced with a turbulence chamber which is similar to those used for simulating conditions of 
observation in astrophysics [25] . These results show that the hybrid encoding is a potential resource for practical long¬ 
distance quantum communication without a shared reference frame, as an Earth-to-satellite communication protocol 
would require. 


Results 


Let us start by considering the effect of turbulence on the unperturbed 0AM eigenstates of a monochromatic 
paraxial field in the cylindrical basis: 

^^(r,0) = (r,0|O = ^(r|Oe*'^ (1) 

where r is the radial coordinate, 0 is the azimuthal angle and I is the orbital angular momentum quantum number. 
For simplicity we will omit the explicit dependence of the wave function on the propagation distance z. We consider 

turbulence in the so-called weak regime, where the random fluctuations of the index of refraction cause only phase 

aberrations in the transverse profile El- This means that an initial 0AM eigenstate transforms as 

= ( 2 ) 

where 0(r, 0) is the random variable associated to the phase perturbation. In general, '0z(r, 0) is a superposition of sev¬ 
eral 0AM eigenstates. Kolmogorov’s theory of turbulence [26] describes the stochastic properties of this perturbation 
by means of the so called coherence function 

= g-6.88x2§(^)i|sin(^)|i^ 

where the subscript “ens” denotes ensamble average. The strength of turbulence in this regime is measured by the 
ratio between the beam waist w and the Fried parameter ro, given by m 


C^z 


ro = 0.185 
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FIG. 1: Basic qubit encoding and decoding procedures considered for hybrid polarization/OAM qubits. A 

linear-polarized single-mode gaussian beam is transformed into a hybrid qubit using a quarter-wave-plate (QWP), half-wave 
plate (HWP) and g-plate device. The decoding procedure works in reverse, and the output beam is sent through a polarizing 
beam splitter (PBS) coupled into a single-mode fiber. 


where A is the wavelength and Cn is the refractive index structure constant. It can be associated with the coherence 
length of the phase perturbations: the smaller the value of ro compared to re, the stronger the effects of turbulence. 

Rotationally invariant photonic qubits are composed of polarization and ^ = ±1 0AM eigenstates. However in 
some cases one might be interested in transmitting higher order hybrid qubits made of polarization and opposite 
0AM values with |/| > 1, as was done in m- The calculations we present here, concerning only the 0AM part, can 
be applied in the future to higher order hybrid qubits. For this reason we consider a two-dimensional 0AM subspace 
spanned by states | — 1) and |/). Moreover, we assume that the detection is such that only these states are registered. 
That is, any perturbation that takes the 0AM outside of the subspace results in a null event, or a loss. This of course 
depends on the detection procedure, and can be true in the case of pure 0AM HO] or hybrid qubits [J. 

Thus, in order to know the properties of the evolved state after turbulence, we have to calculate matrix elements 
of the form {h\'ipi){'ipii\l 2 ) where the 0AM quantum numbers /i( 2 ) and //,(//) take values ±/. The basis elements 
labeled by arable numerals |/i) and I/ 2 ) can be interpreted as the states at the measurement stage. Roman numerals 
denote the initially prepared states |//) and |///) which after turbulence become \ 2 pj) and \^ii) according to equation 
(©• Following the results in na, we proceed to calculate the ensemble average where the non-zero elements can be 
written as 

Cai = {{h\'ipi){‘4’ii\h))ens = J \R{r,z)\^QAi{r,0,d')rdr, (5) 

where R(r, z) is the radial distribution of the unperturbed wave front and the function 0 a/ is the circular harmonic 
transform of the rotational coherence function 

eAi{r,0,0') = J ( 6 ) 

with // — /i — {III — h) = This quantity must be Al = 0, 21 for preparation and detection in the subspace spanned 
by 0AM values Al. 

From equations ^ and ^ one can see that the non-zero elements are 




( 7 ) 
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and 


C2i = \{i^-i\i)\^ = m-i)f- ( 8 ) 

We note that Eq. satisfies the mirror symmetry condition found in [4] for resilience against spatial perturbations. 
Coefficients Co and C 21 depend on re/ro, but we will omit this dependence for notational simplicity. 

Hybrid Qubits. In order to analyze the effect of turbulence on the hybrid qubits, we consider a quantum 
communication scenario without a shared reference frame as in Fig. Logical qubits are encoded into the subspace 
of the polarization/0AM degrees of freedom spanned by the vectors: 


\^)h = \R.l) 
\l)h = \L,-l), 


(9) 


where h stands for “hybrid” and R{L) stands for right (left) circular polarization. One can efficiently produce 
arbitrary superpositions of states in equations ^ by first encoding qubits in the polarization via waveplates and then 
by transferring them to the hybrid space by means of a g-plate, an inhomogeneous birefringent device that couples 
the polarization and 0AM [6l|24]. The g^-plate acts indeed as a universal encoder/decoder, mapping back and forth 
between the polarization qubit |T) = a\R) + P\L) and the hybrid qubit + P\'^)h [I]- 

When I = 1, the hybrid states \0)h and \l)h have null total angular momentum, and are invariant to rotations 
around the propagation axis mi. Hybrid qubits are then sent through a turbulent environment. For the 0AM 
component in equations , after turbulence we have 

I A ^)oam ^ |'0±z)oam5 

where in general {i^i\i^j)oam 7^ 0 even for i which means that orthogonality in the 0AM basis is not preserved 
by turbulence. 

In turn, the effects of turbulence on the polarization degrees of freedom are negligible. So, for hybrid states of light 
in Eq.(|^ after transmission we have 


\0)H^\0^)H = \R)\i^-l) 

li),^liq, = lL)l^,), 


(11) 


where “T” stands for turbulence. Because the polarization is immune to turbulence, orthogonality between logical 
states is preserved even after transmission through turbulence. 

Finally, the detection stage consists of a second g^-plate, wave plates for polarization analysis and a single mode 
fiber. The results 0 and ^ of the last section can be put together to conclude that after the g^-plate, the state 
transforms as [4]: 


a\0^),, + + mW = 0 ) ( 12 ) 

+ higher order I terms, 

where Ph is the probability of the perturbed state to remain inside the logical subspace. By post-selecting the 
fundamental I = 0 spatial mode, the quantum state can be perfectly recovered and encoded back into the polarization 
space. Thus, unity fidelity is in principle possible for the hybrid encoding and any amount of weak turbulence. This 
post-selection can be made in a non-destructive way, by simply sending the photon through a single-mode optical 
fiber. Then, the polarization state can be detected using waveplates and a polarizing beam splitter, as shown in Fig. 
or used for further quantum-information processing. 
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FIG. 2: Experimental setup and calibration of the turbulence chamber a) Experimental setup. A laser beam is 

sent through an encoding stage, consisting of a QWP, HWP and q-plate, used to encode a hybrid qubit. The hybrid qubit 
is then sent through a homemade turbulence chamber, mixing hot and cold air (see text). The decoding stage consists of a 
q-plate and a polarization analysis stage. Light is coupled into single-mode fibers and sent to a photodiode. We also use a CCD 
camera (not shown) to characterize the effects of turbulence. The decoding stage is mounted on a rotation stage, to simulate 
misalignment, b) Characterization of turbulence induced by the chamber given by the ratio of the beam waist w to the Fried 
parameter ro as a function of the power applied to the resistors in the chamber. The insets show snapshots of the beam profile 
taken with a CCD camera with 1ms exposure time. 


Ph constitutes the probability of success of the protocol, and is given by = Cq. Explicitly, using equations (§, 
(§ and we have 


Ph = 


= //i«(r, 


^)l 


2^-6.88x23 (^)3|sin(e-e')|3 


rdrdOdO'. 


(13) 


Experimental Implementation. The experimental setup used to demonstrate the resilience of hybrid qubits 
is sketched in Fig. A linear polarized laser beam with wavelength A = 795 nm propagates first through a 
preparation stage consisting of a quarter wave plate (QWP) a half wave plate (HWP) and a g^-plate where arbitrary 
superpositions of hybrid states can be prepared. A turbulence machine (described in the Methods section) produces 
real air turbulence between the preparation and detection stages. In Fig. [2 b), we observe the relation between 
power at the resistors inside the turbulence machine and the Fried parameter characterizing the turbulence strength. 
Progressive degradation of a gaussian beam is illustrated by 1 ms snapshots for four different Fried parameters. After 
transmission through turbulence, the information is decoded by a g^-plate, followed by QWP, HWP and PBS and the 
light is coupled to a single-mode fiber and detected. 

For 14 different strengths of turbulence, we prepare the 6 states of the set of mutually unbiased bases of dimension 2. 
Fidelity with respect to the initial states is measured at the detection stage obtaining a total of 84 points distributed 
around a mean value of Ph = 0.985 ± 0.006. The experimental results are shown in Fig. Fidelity of transmission 
for the initial states \0)h (orange), \l)h (yellow), |±)/, = {\0)h ± \l)h)/V^ (blue, green), \R)h = (|0)/, - i\l)h)/V^ 
(black), and \L)h = (|0)/i + i|l)/i)/\/2 (red). The effective coupling efficiency decays as the turbulence increases giving 
a measure of the probability of post selection P/^,, shown by the red squares in Fig. The red solid curve is a plot of 
Ph given in equation, ( [l^ , and showing good agreement between experiment and theory. 

To show that the resistance to turbulence is valid in conjunction with the alignment-free character of the hybrid 
qubits, we repeat the experiment for the same set of states but this time we rotate the detection kit along the optical 
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FIG. 3: Experimental results for fidelity and losses a) Success probability (red squares) and fidelity of hybrid states 
(circles) as a function of the turbulence strength w/ro. For the fidelities, each color corresponds to an eigenstate in mutually 
unbiased bases; the grey area encloses the region of dispersion of the experimental points around the mean value 0.985. We can 
observe the resilience to random phase perturbations. b)The rotational invariance of hybrid qubits is preserved in the weak 
turbulence regime. Fidelities do not depend on the angle 0 between sender and receiver’s frames, for ^ = 0.6. The yellow area 
is the region of dispersion of the experimental points. 


axis at several angles. The strength of turbulence was locked at ^ = 0.6. Experimental results are shown in the inset 
of Fig. Again we can conclude the resilience of the qubits by observing a mean fidelity of J^h = 0-98 ± 0.01 taken 
over 30 experimental points, for several different alignment angles. 


Discussion 

We have shown both theoretically and experimentally that rotational-invariant hybrid photonic qubits, encoded 
in the polarization and orbital angular momentum degrees of freedom, are resistant to atmospheric turbulence. The 
robustness is due to both the logical encoding and the decoding procedure, realized with a g^-plate device, which 
provides an intrinsic error-filtering mechanism. Our results show that quantum communication without a shared 
reference frame schemes using hybrid qubits are viable over long distances in free-space. 

We emphasize that our theoretical predictions are based on Kolmogorovs theory of turbulence. Our experimental 
data show excellent agreement with these predictions. The same theory was used in other experiments where thick 
turbulence was simulated by phase screens [28]. Besides, recent experiments in real life situations, as transmission 
of 0AM states of light over a 3 km distance, demonstrate that, even in this cases, Kolmogorov’s theory is still valid 
[23] . This indicates that there are no fundamental obstacles to the application of the present proof-of-principle results 
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to this type of real-life situations. On the other hand, it is known from astronomical observations that stratospheric 
turbulence deviates from the predictions of the Kolmogorov model [29]. Thus, non-Kolmogorov effects will necessarily 
have to be taken into account for Earth- to-satellite quantum communication links. We believe our results constitute 
a basis for future research in that direction. 


Methods 


In order to produce a turbulent medium we construct a machine that causes random fluctuations of the index of 
refraction in the air using a similar design to that used in Refs. HUES]. The light beam is sent through an air 
chamber that is 20 cm long, with total volume 220 cm^. Two cylindrical arrays of resistors with a total electrical 
resistance of 15 Ohms are driven at a maximum power of 290W, heating the air inside the chamber and producing a 
temperature difference of up to AT = 230°C with respect to the room temperature (about 18°C). Two fans behind the 
resistors produce the dynamical turbulence by mixing the air at the room temperature with the heated air, resulting 
in fluctuations of the index of refraction. 

In order to calibrate our machine we place a CCD camera at the waist re of a Gaussian beam after the turbulence 
chamber and observe the average broadening of the Gaussian profile after propagation in the turbulence chamber, as 
compared to the unperturbed beam. For a given value of the power applied to the resistors, we record a 20 s video 
using the GCD camera. Summing over all frames of the video, we measure the width of the Gaussian distributions 
with and without turbulence. Following Refs. [niiiQ!, the ratio w/vq can be obtained via the equation 


w 

ro 



(14) 


where Wt is the width of the beam after turbulence. This allows us to measure values from 0 to a maximum value of 
— = 1.4 ± 0.2, corresponding to a wide range of the turbulence regime. 
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